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Abstract 

We study the short-time dynamics of systems that develop "quasi long-range 
order" after a quench to the Kosterlitz-Thouless phase. With the working hy- 
pothesis that the "universal short-time behavior", previously found in Ising- 
like systems, also occurs in the Kosterlitz-Thouless phase, we explore the scal- 
ing behavior of thermodynamic variables during the relaxational process fol- 
lowing the quench. As a concrete example, we investigate the two-dimensional 
6-state clock model by Monte Carlo simulation. The exponents governing 
the magnetization, the second moment, and the autocorrelation function are 
calculated. From them, by means of scaling relations, estimates for the equi- 
librium exponents z and n are derived. In particular, our estimates for the 
temperature-dependent anomalous dimension n that governs the static corre- 
lation function are consistent with existing analytical and numerical results 
and, thus, confirm our working hypothesis. 

PACS: 64.60.Ht,75.40.Gb,67.40.Fd,02.70.Lq 
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I. INTRODUCTION 



As discovered by Janssen et al. 0, universal short time behavior (USTB) in critical 
dynamics occurs after quenching a spin system with nonconserved order parameter (model 
A H|) from a high-temperature initial state to the its critical point. In its most pronounced 
form the USTB manifests itself in the time dependence of the order-parameter expectation 
value. Starting from a small average initial magnetization m , the order starts to grow during 
a nonequilibrium (though universal) period shortly after the quench, before the decay to the 
equilibrium value zero sets in. The initial increase is governed by a power law, 

m ~ mo t e , (1.1a) 

and the short-time exponent 9 is given by [||,[| 

6>=*LZ^, (1.1b) 

i.e., it is determined by the difference between the anomalous dimensions t]q and t] of the 
initial order-parameter field and the equilibrium magnetization, respectively, (z denotes the 
dynamic equilibrium exponent.) 9 (or, equivalently 770) is a new exponent, not expressible 
in terms of known static and/or dynamic exponents PJ. Thus, the initial increase will only 
occur under the condition that 770 > rj, which, for instance, is fulfilled in the case of the Ising 
model for spatial dimensions 1 < d < 4. 

Recently significant progress has been made in quantitative studies of the USTB by 
means of Monte Carlo (MC) simulations. It turned out to be possible to observe the short- 
time power law ( |l.laj) in surprisingly small systems and after relatively short simulation 



times As the exponents governing the short-time behavior of thermal averages can 

be expressed in terms of scaling relations among short-time and equilibrium exponents, it 
is possible to extract even static exponents from a relatively early stage of the dynamics, 
whereby the problem of critical slowing down is circumvented to some extent 

As pointed out by Janssen |§, the qualitative mechanism underlying the USTB is the 
following: After the quench the initially small correlation length starts to grow towards its 
equilibrium value. During the early stages, when fluctuations are still short-ranged, a mean- 
field description of the process should hold. Since in Ising-like systems T C MF is always larger 
than the real T c , it is plausible that shortly after the quench a spontaneous magnetization 
will evolve, pointing in the direction of the initial magnetization m . At some point, when 
£(t) has grown to a certain size, this increase will come to a halt, and afterwards the decay 
to the equilibrium value m(t — > 00) = will begin. 

So far in both analytical and numerical calculations the USTB has only been studied in 
Ising-like systems, with a conventional critical (or tricritical) point fl],3,[5HTT|l • However, the 
phenomenon (and the applications built thereon) should not only occur at a critical point, 
but, more generally, in situations where scale-invariant long-ranged fluctuations evolve after 
a quench from a high-temperature initial state (or after any other comparable change of the 
state of the system). In the following this scenario is tested for a quench to the Kosterlitz- 
Thouless (KT) phase? In two-dimensional systems with 0(2) symmetry (or variants of it) 
below the KT transition temperature Tkt, one finds a range of temperatures with a power- 



law decay of the static correlation function |12|JT3|| . Examples are the XY model and the 
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p-state clock model for p > 5. Additionally, mean-field theory predicts for these systems 
(erroneously) a conventional critical point, and one finds generally T C MF > Tkt- Thus, taking 
into account the foregoing discussion, the USTB should also be observable after a quench 
to the KT phase, and it is interesting to see whether the methods described in Refs. |5|-|8| 
can be applied, i.e., whether for instance the equilibrium exponent rj can be extracted from 
the short-time relaxational behavior. 

The rest of this paper is organized as follows. In the next Section we recall analytical 
and numerical results for the clock model. In Sect. Ill we derive the scaling behavior of the 
relevant thermodynamic quantities, taking into account the anomalous scaling of the initial 
magnetization. Our own simulations have been carried out for the p = 6 clock model. The 
numerical results are presented in Sect. IV. 



II. THE CLOCK MODEL 



According to the Mermin- Wagner theorem ||14|| , the two-dimensional XY model cannot 



undergo a conventional para-/ferromagnetic phase transition. However, as shown by Koster- 



litz and Thouless [12]], there exists a transition between a state with unbounded and bounded 
vortices at T = T K t- 

A model that in many respect interpolates between the Ising model and the XY model 
is the p-state clock model, defined by the Hamiltonian 

H = - J $i ■ Sj = - J cos - tij) , i?j = — Hi , (2.1) 



<ij> <ij> 



P 



where e {0, 1, ... ,p — 1} and < . > denotes the sum over nearest neighbors. This spin 
system is equivalent to the XY model for p — > oo and to the Ising model for p = 2. It was 
first shown by Jose et al. []T3 1 and later on verified by numerous MC studies [pT5] pOH that the 



p-state clock model has a single Ising-like critical point for p < 4. For p > 5, on the other 
hand, there exist two transition points, with temperatures T\ and T 2 (> 7\), respectively. 
For T > T2 a paramagnetic high temperature phase exists, and for T < T\ the system 
becomes ferromagnetic. For T\ < T < T2, however, a KT phase with "quasi long-range 
order" emerges, where the correlation function behaves as C(r) ~ |r| _T? , with the anomalous 
dimension rj = f](T) continuously depending on the temperature JTJ]. Both transitions are 
of KT type, i.e., the correlation length diverges exponentially as 7\ (T 2 ) is approached from 
below (above) [|T3|1 . 

Our MC simulations — the results are presented in Sect. IV — were done for the p = 6- 
state clock model. The case p = 6 was studied in great detail in the literature, partly 
motivated by the fact that in real two-dimensional systems one expects anisotropies with 
six-fold symmetry due to the lattice structure [O. Approximate analytic results for the 



transition temperatures and the exponent Tj(T) have been obtained by Jose et al. [13 
While T2 was found to be identical with the KT transition temperature of the XY model, 
T/J=0.95, T\ was predicted to vary as T\j J = 47r 2 /1.7p 2 (which means T\j J = 0.645 for 
p — 6). (Here and in the following we set the Boltzmann constant ks = 1.) The exponent 
77 was found to vary between r](T 2 ) = 1/4 and i](Ti) = 4/p 2 (which means rj(Ti) = 0.111 for 
p = 6). 
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Among the MC analyses mentioned above, we used the results of Challa and Landau 
for comparison with our own data. They found for the transition temperatures T 2 /J = 
0.92(1) and T%/ J = 0.68(2) and for the corresponding exponents rj(T 2 ) = 0.275 and r)(Ti) = 
0.100. (No errors for rj(T) are given in Ref. ||18|| .) Both the results for the transition 
temperatures and the exponents are consistent with the analytic results of Ref. |L3] . 

Dynamic MC simulations for the 6-state clock model were performed by Kaski et al. 
[T7J , but in this work the relaxational behavior in the ferromagnetic phase was studied, with 
emphasis on the growth of the domain size after a quench. Thus, concerning the temperature 
ranges, this work is complementary to ours. 

In order to comply with the requirements for the occurrence of the USTB discussed in 
the Introduction, the system has to have a mean-field T C MF above the temperature to which 
the quench is carried out. In the case of the clock model it is straightforward to show that 
there exists a T C MF for all p > 2 and that T C MF / J = 2d for p = 2 (Ising case) and T C MF / J = d 
for p > 2, where d denotes the spatial dimension. Comparison with exact or numerical 
results for the clock models shows that the upper transition temperature (which is a critical 
point for p — 2, 3, 4) is below T C MF for all p > 2. E.g. for p = 6, the upper transition 
temperature T 2 /J ~ 0.92 is well below T C MF / J = 2. 



III. SCALING ANALYSIS 

In the KT phase, i.e. for all temperatures below Tkt, the bulk equilibrium correlation 
length is formally infinite. After a sudden quench from a high-temperature initial state 
to the KT phase, the initially small correlation length £(t) starts to grow, and the large- 
scale fluctuations gradually develop. Let us first discuss relaxational processes in Ising-like 
systems. In this case, after a short period of nonuniversal behavior following the quench, 
the growth of the correlation length is described by the universal power law, 

m-t l,z , (3.i) 

and in the bulk system £(t) becomes "macroscopic" for large times. Moreover, in the uni- 
versal regime £(i) is the only macroscopic time scale and allows to characterize the time 
dependence of correlation functions completely. In the finite systems, ( |3.1| ) holds as long as 
£(t) < L (where L is the linear system dimension). In the following this scenario will be 
called simple (dynamic) scaling. 

In the KT phase vortices are generated in the course of the quench and potentially 
disturb the simple scaling picture |2~T|-F2"3|| . Deviations from the simple scaling picture are 



expected mainly for intermediate times. As we shall see from the evaluation of the MC data, 
modifications compared to the Ising dynamics really occur in certain quantities in form of a 
relatively slow approach to the power-law form. Below we review the simple scaling scenario 
concerning the quantities studied in our MC simulation, with emphasis on the (expected) 
consequences of the USTB. Modifications due to vortices are discussed in some more detail 
in the context of the MC simulation in Sect. IV. 
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A. Bulk behavior 



Consider first the expectation value of the magnetization density |24]. Suppose in the 



high-temperature initial state an external magnetic field has generated a small (initial) 
magnetization m , pointing in the i? = direction of the xy plane. Then, after the quench, 
the expectation value of the order parameter, defined by 

m(t) := L- 2 (J2cos$i) , (3.2) 

i 

should remain in the direction of the initial magnetization. In the scaling regime it should 
further satisfy the homogenity property PS] 

m(t, mo) « l^ 2 m{l z t, r" o/2 m ) , (3.3) 

where / is a rescaling parameter, and the exponents rj, r/o were introduced in Eqn. QLlaQ . 



From Eqn. ( |3~3|) it follows that, asymptotically, the magnetization behaves as 

m{t, mo) ~ r^ 2 * f(m * W2 *) , (3.4) 

where T{x) is a universal scaling function. As an immediate consequence of 7/0 7^ 77 in Eqn. 
Q3.3|) , the rescaled (or renormalized) mo is no longer the initial value of m(t). However, one 
still has asymptotically m oc m for t — > Thus the scaling function J^(x) in ( |3.4|) 

behaves ~ x for a; — > 0, and, as a result, the short-time behavior is described by ( |1 . 1|) . 

Further, from the combination in the argument of T in (|3.4j ) one realizes that mo also 
gives rise to a novel time scale, t ~ m 2z ^° . While for t < t the USTB ( |l.la| ) is observed, 
the system crosses over to the nonlinear decay m ~ t~ v ^ 2z for t > t , which implies that 
J-{x) — > const, for x — > 00 and that memory of the initial value mo is lost for t > t - 

Can we expect to observe the USTB even for m = 0? Obviously, in this case the 
magnetization vanishes identically for t > 0. However, as demonstrated in Ref. ||, the 
USTB leaves its fingerprint on the relaxational behavior of the autocorrelation function 

A{t) = (s t (t) ■ si(0)) , (3.5) 

which measures the change the system has undergone during the time t after the quench. 
Without going into the technical details, we here just present Janssen's result @: 

A[t) ~ r x/z , with A = 2 - z6 , (3.6) 

i.e., the short-time exponent 9 also appears in A(t) and leads to a deviation from the naive 
expectation A = d, surprisingly for all £ > 0. It was this quantity where the USTB actually 
was first observed in a MC simulation pSfl . 

The third quantity that will play a role in the following is the second moment (or variance) 
of the magnetization. It is defined by 



m 



(2) 



(t) = L- 2 (j:m-m)- w 



■1.1 



A straightforward analysis reveals that m^ 2 '(t), again for vanishing initial magnetization, 
takes the scaling form 
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m (2) (t)~t ft , with b=^—^-, (3.8) 



i.e., the USTB does not influence this quantity P? |. 

As suggested by Schiilke and Zheng ||, the three exponents 9, A, and b, governing the 
bulk behavior of m, A, and mS 2 \ respectively, suffice to determine the equilibrium quantities 
7] and z from the USTB. The starting point is the short-time exponent itself, which can 
be obtained from simulations with small, nonvanishing itiq. Then, from simulations with 
mo = 0, the exponents A and b can be estimated simultaneously. Utilizing the scaling 
relation ( [3.6D , z can be determined, and, in turn, with relation ( |3.8|) eventually the scaling 
dimension r]. 

Our main hypothesis for the present work is that also in the KT phase, analogously to the 
situation in the critical Ising model, the operator dimension of the initial field, 7/0, is different 
from r], and that, as a consequence, the scaling picture developed above also describes the 
short-time dynamics in the KT phase. This is neither obvious nor by any means rigorously 
proved. At the present stage it is merely a working hypothesis, which below will be tested by 
means of MC simulations. Especially for quenches to the KT phase, we expect all exponents 
to depend on the respective final temperature to which the quench is carried out. 



B. Finite-size scaling 

Starting from Eqn. (|3.3|) , all the results about the scaling behavior of the thermodynamic 
quantities described so far hold in the bulk system. If we are dealing with a system of finite 
size — and certainly this is the case in the MC simulation — the power laws ( |l.la|) , ( p.6|) 



and ( |3.8| ) hold as long as the time-dependent correlation length £(t) ~ t x l z is smaller than 
the system size L. A finite-size scaling analysis reveals how the thermodynamic quantities 
behave during the later stages of the relaxational process P,p5[. 

Taking into account the finite system size, the analog to (|3.3|) reads 

m(t, m , L) « V l2 m{lH, r w/2 m , XL) . (3.9) 

Removing the rescaling parameter / by setting / ~ L, one obtains 

m(t, m , L) ~ L~' n/2 G(t/L z , m L no ^ 2 ) , (3.10) 

where Q{x,y) is a universal scaling function. As seen from ( |3.10|) , the size L of the system 
gives rise to another characteristic time scale ti ~ L z (besides the initial time scale to), 
the well-known finite-size scale f29| . In the KT phase it turns out temperature dependent. 
Further, it was shown in Ref. |§ that Q(x, y) ~ A(y)e~ x for x — » 00, i.e., the magnetization 
decays linearly when t 3> L z . The scaling function A(y) in ( |3.10|) describes the universal 
dependence of late stages of the relaxational process on the initial magnetization m . This 



dependence is special to the finite-size system |28|.In the bulk memory of the initial condition 
is lost for t — > 00. 

For arbitrary x and y — > 0, the function Q(x,y) behaves as Q ~ yB(x), with a scaling 
function B. In this limit the maximum of the magnetization profile has the scaling form 

m max ~ LM/ 2 m B(t ma jL z ) . (3.11) 
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Thus, the finite-size scaling behavior of m max allows in principle to determine directly the 
difference rj — rj. 

Also for the other quantities discussed in Sect. III. a the asymptotic behavior is modified 
when £(t) > L. Firstly, the autocorrelation function ( |3.5|) also decays linearly for t > t^. 
Secondly, the second moment takes the scaling form 

mW(t,L)~t b J(t/t L ). (3.12) 

In order to reproduce ( |3.8|) , the scaling function J{x) has to approach a constant for x — > 0. 
For x — > oo, J(x) ~ x~ b , such that in the long-time limit the equilibrium result rnP^ ~ L 2 ~ v 
is obtained. 

Again, in the KT phase we expect to observe the finite-size scaling scenario presented 
above, with the characteristic scale ti depending on the temperature. 



IV. MONTE CARLO ANALYSES 
A. The method 



We used the usual time-dependent interpretation of the sampling algorithm |3(J, with 
sequential updating and nonconserved total spin. We employed the heat-bath algorithm, as 
described in detail for the clock model in Ref. p0|. It is equivalent to Glauber dynamics 



30fl . All simulations were carried out on square lattices with varying linear dimension L 



with periodic boundary conditions. 

As the initial configurations we used equilibrium states for uncoupled spins in an external 
magnetic field h = 2mo (with mo = for the autocorrelation function and the second 
moment). Time-dependent thermal averages were approximated by summing over a number 
of N histories, N varying between 1000 and 10000 depending on the size L. In order to 
obtain a measure for the error, the N histories were divided in 10-20 runs, from which 
statistically independent exponent "measurements", their mean value, and the statistical 
error of the mean value were obtained. 

As mentioned in the Introduction, an attractive feature of the USTB in Ising-like sys- 
tems is that relatively relatively small lattices and short simulation times turned out to be 
sufficient extract reliable results for exponents ||[7|,||] . However, the estimates obtained will 
in general to some extent depend upon the lattice size L and, in the case of 9, upon the 
initial magnetization mo- In the case of the clock model, it turned out that one has to go to 
somewhat larger lattices and, especially in the case of the autocorrelation function, the time 
until the power-law behavior is assumed is increased compared to the Ising or Potts model. 



B. The short-time exponent 6 and the linear relaxation time 

Firstly we discuss our data for the magnetization and the determination of the short-time 
exponent 6 from them. In order to get an impression of the time dependence of the order 
parameter, we have calculated relaxation profiles of the 6-state clock model for L = 30 and 
m = 0.05 in the temperature range 0.5 <T/J< 1.1, i.e., the lower limit lies well below Ti 
in the ferromagnetic phase and the upper limit above T 2 in the paramagnetic phase. 
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Fig 1.: Magnetization of the 6-state 
clock model as a function of time 
(in units of Monte Carlo steps per 
spin) in semi-logarithmic represen- 
tation for L = 30 and m,Q = 0.05 for 
various temperatures in, below, and 
above the KT phase. From bottom 
to top the corresponding tempera- 
tures are T/J = 1.1,1.04, 1.0, 0.95, 
0.92, 0.84, 0.76, 0.72, 0.68, 0.64, 0.6, 
0.5. 



t (MCS) 




Fig 2.: The same profiles as in 
Fig. 1 in double-logarithmic repre- 
sentation, emphasizing the initial 
power-law growth of the magnetiza- 
tion. 



t (MCS) 



The results are displayed in semi- logarithmic form in Fig. 1. and in double-logarithmic 
form in Fig. 2. For temperatures T/J < 1.1 the magnetization has the tendency to increase 
initially. For the highest temperature represented in the figure, this increase lasts only for 
about 10-20 updates of the lattice, but it becomes more and more pronounced at lower 
temperatures. As one can see from Fig. 2, a clearly identifiable initial power-law regime 
at about 10 < t ;$ 100 can be observed for T/J < 0.95. (The time is always given in 
Monte Carlo steps per spin.) All curves in the temperature range 0.64 < T/J < 1.1 have a 
maximum and thereafter decay exponentially or "linearly" ~ e~'/' L . 

It is instructive to extract the linear relaxation time from the data in Fig. 1. It depends 
on both the system size and the temperature; results for ^l(T) for L = 30 are displayed in 
Fig. 3 (squares for p — 6). The temperature range of the KT phase as obtained by Challa 
and Landau fl8| is indicated in the figure. Obviously, in this quantity the signature of 
the KT phase is a plateau-like structure, roughly located between Ti and T 2 of Ref. [JTS . 
Coming from the high-temperature side, t L increases up to about T/J = 0.95. Then t L 
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remains almost unchanged down to T/J = 0.75 and starts to grow rapidly when going to 
lower temperatures. For T/J < 0.6 the magnetization more or less approaches a constant 
within the observation time t < 1200, which means that becomes effectively infinite. The 
latter is due to an exponentially growing tunneling time in the ferromagnetic phase [j30|j . 
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Fig 3.: The linear relaxation times 
ti for p = 6 (squares) and p = 4 
(circles) function of tempera- 

ture for L = 30, mo = 0.05. Both 
the (exact) critical temperature of 
the 4-state clock model, T c /J = 
1.13.., and the upper and lower tran- 
sition temperatures T\/J = 0.68 
and T2/J = 0.92 determined Challa 
and Landau Jyl are indicated in the 
figure. The error bars are approxi- 
mately of the same size as the sym- 
bols. 



For comparison we calculated a number of curves for the 4-state clock model in the 
vicinity of its (exact) critical temperature T c /J = 1.1345.... The corresponding linear 
relaxation time is also shown in Fig. 3 (circles). It is steadily increasing for decreasing 
temperature, as one would have expected in the vicinity of a critical point, and no plateau 
is seen. 

In order to determine the short-time exponent 9, we calculated magnetization profiles for 
the lattice size L = 300. As also found for other models ||, the outcome for the short-time 
exponent 9 varies significantly with m^. From the renormalization-group analysis it is clear 
|]J3|] that the pure scaling form ( |l.la| ) is assumed for m — > only. Unfortunately, it is hard 
to obtain reliable results for 9 in this limit, as for mo = the signal vanishes completely, 
and for decreasing mo it becomes more and more blotted out by the noise. Hence, as in Ref. 
||, we calculated profiles for a number of initial values mo and then determined 9 by fitting 
a power law to the short-time regime. As the fit interval we chose [10, 300]. The results for 
9 as a function of mo for the temperatures T/J = 0.6, 0.68, 0.76, 0.84, 0.92 are shown in 
Fig. 5. The values of the short-time exponent obtained by extrapolating linear fits to m = 
are listed in the Table. 

Some representative profiles for L = 300 and m = 0.01 are presented in Fig. 4 (solid 
lines). They show that in the given temperature range the initial increase is consistent with 
a power law. As seen from the other curves displayed in Fig. 4, this is drastically different 
at low temperatures. E.g. for T/J = 0.2 (upper dashed curve), deep in the ferromagnetic 
regime, the time dependence of the magnetization is more complicated, and only for t ^ 200 
power-law behavior is assumed. This is weakly indicated already at T j J = 0.6 (lower dashed 
line). 
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Fig 4.: Initial behavior of the order 
parameter for L = 300 and mo = 
0.01 for (from bottom to top) T/J = 
0.92, 0.84, 0.76, 0.68, 0.6, 0.2. The 
solid curves represent profiles within 
the KT phase. 




Fig 5.: Results for the short-time 
exponent 9 as a function of uiq 
for four temperatures T/J = 0.68 
(squares), 0.76 (full circles), 0.84 
(triangles pointing up), and 0.92 
(triangles pointing down). The re- 
spective linear fits, from which the 
extrapolations to mo = were ob- 
tained, are also displayed. 
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C. The autocorrelation function and the exponent A 



The autocorrelation function A(t) as defined in Eqn. ( |3.5| ) is governed by the exponent 
X/z. Exemplary results for A(t) in the KT phase for L = 300 (and mo = 0) are displayed 
in Fig. 6 (solid lines). Even from this representation one can tell that the curves are slightly 
curved up to t ~ 100 — 200 and that, if at all, they only relatively slowly approach a power 
law. This is the reason why the exponent X/z in (|3.6| ) cannot be determined on smaller 
lattices with, say, L = 100. In this case A(t) crosses over to the finite-size (linear) decay 
before a reliable estimate for X/z can be extracted. For comparison we have also plotted 
the autocorrelation function for the critical Ising model in Fig. 6 (dashed curve). 

In order to overcome this obstacle, we calculated the autocorrelation function for L = 300 
up to t = 1000. Then X/z was determined by fitting the data in an interval [ti,t 2 ] within 
the power-law regime. To determine the minimal t\ (called delay time in the following), 



10 



we calculated the difference AA between the data and the fit for a given interval. While 
the upper limit t<i ~ 1000 was fixed, the lower limit t\ was increased until (besides random 
fluctuations) AA did not show any systematic deviation from the fit within [ti,/^]- 




Fig 6.: The autocorrelation func- 
tion defined in ( |3.5| ) for L = 300, 
m = for T/J =0.92, 0.84, 0.76, 
0.68, 0.6, 0.2 (upper six curves from 
bottom to top) The dashed curve 
represents the data for the critical 
Ising model. 
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Consider first the temperature T/J = 0.68. Numerical results for AA are shown in Fig. 7 
(two upper curves). For the uppermost curve the fit interval was [200, 1000]. It is obvious 
that, besides the random short-scale fluctuations, the data systematically deviate from the 
fitted power law; for 200 < t < 300 the data are below and for 400 < t < 600 they are 
systematically above zero. In contrast, a fit to the interval [500, 1000] (second curve from 
above in Fig. 7) shows no systematic deviation, and, thus, for this temperature the power- 
law regime seems to be reached for t ~ 500. 
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Fig 7.: The difference AA(t) be- 
tween the data for the autocorre- 
lation function and the power-law 
fit obtained for T/J = 0.68 and fit 
intervals [200, 1000] (upper curve) 
and [500,1000] (middle curve). The 
bottom curve shows AA for the 
critical Ising model and fit interval 
[80,1000]. 
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In Fig. 7 we have also plotted AA(t) for the critical Ising model with fit interval [80, 1000]. 
Although the delay time t% ~ 80 is much shorter than for the clock model, it is still relatively 
large compared with m(t) or m^ 2 '(t), where one finds power-law behavior for t > 10. 
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Fig 8.: AA(t) for T/J = 0.92, 
0.84, 0.76, 0.6, and 0.2; the re- 
spective fit intervals are [200,1000], 
[450,1000], [650,1000], [400,1000], 
and [500,1000]. 
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In Fig. 8 the difference AA(t) is displayed for various temperatures within and below the 
KT phase. The estimated delay times are 500, 400, 500, 650, 450, and 200 for T/J = 0.2, 
0.6, 0.68, 0.76, 0.84, and 0.92, respectively. Thus, t x is relatively short at T 2 , increases 
when T is decreasing, and appears to become smaller again for T < T\. However, presently 
our data are not accurate enough to allow for definite statements about the temperature 
dependence of t\. 

The results for X/z are shown in the Table. The slow approach the power- law behavior 
of A(t) is a new feature not seen in the Ising model. Presumably it is due to vortices which 
are abundantly generated during the quench and slowly annihilate during the ensuing relax- 
ation. The influence of vortices on the dynamics of the XY model was discussed by several 
authors in the literature. Vortices in general give rise to new length scales f2l"lj23|j , and, 
as a consequence, the simple scaling picture discussed in Sect. II breaks down. It was also 
shown that in certain quantities a slow (logarithmic) approach to the power-law behavior 
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occurs [21] , 22] , especially in those which are substantially influenced by short wavelengths 
fluctuations of the order-parameter field. However, with the help of the existing literature 
it was not possible for us to get a complete understanding of the observed dynamics at 
intermediate time scales. Because we were mainly interested in the asymptotic power laws 
and the exponents governing them, we did not pursue this point any further. 



T/J 


e 


X/z 


b 


z 


V 




0.92 


0.254(5) 


0.670(10) 


0.816(12) 


2.16(4) 


0.23(6) 


1.10(3) 


0.84 


0.273(6) 


0.644(19) 


0.832(10) 


2.18(6) 


0.19(7) 


1.19(4) 


0.76 


0.297(2) 


0.622(11) 


0.842(12) 


2.18(3) 


0.17(5) 


1.29(2) 


0.68 


0.314(2) 


0.577(6) 


0.856(11) 


2.24(2) 


0.08(4) 


1.41(2) 


0.60 


0.355(5) 


0.517(6) 


0.838(11) 


2.29(3) 


0.08(5) 


1.63(3) 


Ising 


0.191(2) 


0.730(6) 


0.806(2) 


2.172(6) 


0.25 


0.83(1) 



Table 1: Monte Carlo estimates for the exponents of the two-dimensional 6-state clock model for 
temperatures within and slightly below the KT phase. For comparison the respective exponents of 
the critical Ising model (at T c /J = 2.2691 . . .) are also displayed. In the latter case 9 and z were 
taken from Ref. ||, then X/z, b, and ijo — r) were obtained with the help of Q, Q, (|LTb| ) using 
the exact r\ = 0.25. 



D. The second moment and results for z and rj 

The short-time behavior of the second moment mS 2 \t) defined in ( |3.7|) is much less 
problematic than that of A(t). Some representative curves are shown in Fig. 9. As before, 
the results for the KT phase are the solid curves. The power-law behavior is found for 
t > 20, and the exponent b in ( |3.8j ) was determined in the interval [20, 500]. The results are 
listed in the Table. 

Obviously, b and because of (|3~8| ) the combination (2 — rj)/z do not depend much on 
the temperature. Taking into account the errors, the outcome is even consistent with a 
temperature-independent exponent, although a weak decrease for increasing temperature 
seems also possible. 

As can be further seen from Fig. 9 also for the critical Ising model (dashed curve) this 
exponent is not much different. From our own data displayed in Fig 9 we determined b = 
0.81(1). Taking alternatively the exact rj = 0.25 and the literature value z = 2.172(6) ||, 
one finds with the help of ( |3.8D b = 0.806(2). Finally, we use the scaling relation ( ^.6| ) and 
( P-8| ) to get estimates for the dynamic exponent z(T), the anomalous dimension r](T), and 
the difference 7] — rj. The results are also displayed in the Table. 



13 



1000 



m (t) 



100 




Fig 9.: The second moment de- 
fined in for L = 300, m = 
and T = 0.92, 0.84, 0.76, 0.68, 0.6 
(lower five curves from bottom to 
top). The dashed curve shows the 
data for the critical Ising model. 
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V. SUMMARY AND DISCUSSION 

We investigated the short-time relaxational behavior of the two-dimensional clock model 
after a quench from a high-temperature initial state to the KT phase. Motivated by the 
situation in Ising-like systems, we assumed that an independent initial exponent also exists 
for the clock model, that it determines the short-time behavior in the KT phase, and that the 
static exponent rj(T) can be determined from the initial stage of relaxation processes. Under 
the assumption that the dynamics is asymptotically governed by only one (macroscopic) 
length scale, the time-dependent correlation length £(t) ~ t 1 /*, we explored the consequences 
of the anomalous short-time behavior. Numerical results for the magnetization, its second 
moment and the autocorrelation function were obtained by Monte Carlo simulation. 

A new feature that is not found in the Ising model is a relatively slow approach to a 
power law of the autocorrelation function. Presumably, this is due to vortices and additional 
length scales associated with them. From the existing literature we were not able to get a 
full quantitative understanding of this phenomenon, in particular in view of the relatively 
fast approach to the power law observed in the magnetization and its second moment. 

Another observation concerns the exponent b of the second moment. It turns out to be 
either weakly depending on the temperature or even constant. It would be interesting to 
study the temperature dependence of this quantity also for other values of p and for the XY 
model (p — > oo). 

Our Monte Carlo estimates for rj(T) listed in the Table are consistent with earlier ana- 
lytical and numerical work [|l|,[]l|]. Since r\ is a derived quantity, its error is relatively large, 
and with the present procedure it would be hard to obtain a much smaller value. However, 
we believe that these results confirm our basic hypothesis that the universal short-time be- 
havior, previously known from Ising-like systems, also occurs in the KT phase. In particular 
this implies that there exists a new exponent r)o(T), the anomalous dimension of the initial 
magnetization, that in general is different from the exponent rj(T) of the static correlation 
function. This difference gives rise to the characteristic short-time behavior of the magne- 
tization (|l.la|) and also leaves its fingerprint on the time dependence of the autocorrelation 
function A(t) as expressed in (|3.6|). 
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